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Nematic  liquid  crystal  polymers  and  nanocomposites  have  wide-ranging  applications  in  modern 
technology  including  display  devices,  ultra-fast  switches,  high  strength  fibers,  and  materials  with 
enhanced  multi-functional  properties  including  thermal,  dielectric,  electrical,  and  barrier  properties. 

In  this  review  article  we  provide  an  overview  of  selected  mathematical  issues  and  numerical  simu¬ 
lations  of  nematic  liquid  crystal  polymers  and  rigid  rod  nanocomposites.  Some  open  questions  will 
be  addressed. 
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1.  INTRODUCTION 

Macromolecular  materials,  in  particular,  liquid  crystalline 
materials,  have  a  wide  range  of  applications  in  ultra-fast 
video  displays,  temperature  and  pressure  sensors,  chro¬ 
matography  and  smart  fluids  for  brakes  and  clutches,  high- 
strength  fibers  and  many  more.95  One  prominent  example 
is  that  nematic  liquid  crystal  polymers  have  been  spun  into 
high  strength  fibers  and  then  manufactured  into  special 
airbags  that  cushioned  the  successful  landing  of  NASA’s 
highly  publicized  missions  to  Mars.22 

Peter  Palffy-Muhoray  has  given  many  interesting  exam¬ 
ples  of  liquid  crystals  in  our  life  in  Physics  Today :88  “Liq¬ 
uid  crystals  are  all  around  us:  in  high  strength  plastics, 
snail  slime,  laundary  detergents,  textile  fibers  such  as  silk 
and  Kevlar,  crude  oil,  insect  wings,  mineral  slurries,  lip¬ 
stick,  Bose  Einstein  condensates,  and  the  mantles  of  neu¬ 
tron  stars.  We  eat  them  as  aligned  molecules  in  gluten  and 
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drink  them  as  phospholipids  in  milk  where  they  stabilize 
fat  globules.  In  our  bodies  they  transport  fats,  make  up 
cell  membranes,  and  affect  the  functioning  of  hair  cells  in 
the  inner  ear,  and  even  of  DNA.”  So,  what  are  liquid  crys¬ 
tals?  Liquid  crystallinity  refers  to  a  state  of  matter  that  is 
intermediate  between  the  crystalline  solid  and  the  amor¬ 
phous  liquid.  As  a  rule,  a  substance  in  a  liquid  crystal  state 
is  strongly  anisotropic  in  some  of  its  properties  and  yet 
exhibits  a  certain  degree  of  fluidity.  When  the  concentra¬ 
tion  of  the  rodlike  molecules  is  high  (lyotropic  systems) 
or  when  the  temperature  is  low  enough  (thermotropic  sys¬ 
tems),  the  molecules  tend  to  be  aligned  with  each  other  and 
the  system  is  in  a  nematic  phase.  Otherwise,  the  molecule 
orientations  are  randomly  distributed  and  the  system  is  in 
an  isotropic  phase.  Figure  1  illustrates  the  arrangement  of 
rod-like  molecules  in  both  the  isotropic  (or  random)  phase 
and  the  nematic  phase.  Actually  there  are  many  anisotropic 
mesomorphic  phases  or  mesophases  (meaning  intermedi¬ 
ate  phases)  with  the  nematic  phase  being  the  simplest  one. 
The  orientational  order  presented  in  liquid  crystals  is  usu¬ 
ally  quantified  by  an  order  parameter  while  the  orienta¬ 
tion  of  a  rodlike  molecule  is  specified  by  a  unit  vector 
along  its  axis.  The  nematic  liquid  crystals  have  a  high 
degree  of  long-range  orientational  order  of  the  molecules 
but  no  long-range  positional  order.  This  solid-liquid  dual¬ 
ity  leads  to  many  anomalous  behaviors  of  liquid  crys¬ 
talline  materials.  Liquid  crystalline  polymers  (LCPs)  are 
materials  with  long  and  stiff  macromolecules,  whose  con¬ 
formation  and  flexibility  interact  with  fluid  flows,  applied 
electric  or  magnetic  fields,  and  thermal  conditions.  There 
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are  two  types  of  such  materials,  one  called  lyotropic  which 
is  a  solution  that  undergoes  a  phase  transition  between 
isotropic  and  nematic  states  when  the  concentration  is 
increased  sufficiently,  and  the  other  called  thermotropic 
which  are  melts  that  exhibit  the  isotropic  to  nematic  phase 
transition  when  the  temperature  is  decreased  to  a  critical 
temperature.  The  distinction  between  lyotropic  and  ther¬ 
motropic  systems  is  not  crucial  for  mathematical  studies. 
Products  with  high  tensile  modulus  are  manufactured  from 
these  macromolecular  materials  due  to  preferred  internal 
molecular  orientational  order.  Understanding  the  flow  of 
LCPs  and  the  effects  of  external  fields,  in  particular,  the 


flow-orientational  order  coupling  during  material  process¬ 
ing,  is  crucial  to  predicting  new  materials  with  desired 
properties  and  to  quality  control  in  LCP  material  process¬ 
ing.  Excellent  books  on  this  topic  include  Bird  et  al.,3,4 
Doi  and  Edwards,23  de  Gennes,20  Larson,67  and  Donald 
et  al.22  A  detailed  overview  on  liquid  crystals  can  be  found 
in  Ref.  [96], 

The  initial  theoretical  basis  for  the  rigid  rod  liquid  crys¬ 
talline  polymer  is  attributed  to  two  Nobel  Prize  laureates — 
Lars  Onsager  and  Paul  J.  Flory.  Both  the  Onsager  and 
Flory  theories  are  statistical  theories  based  on  the  assump¬ 
tion  that  liquid  crystal  phases  arise  from  purely  steric 
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(a) 


Fig.  1.  The  arrangement  of  rod-like  molecules  (a)  in  the  isotropic  phase; 
(b)  in  the  nematic  phase. 

interactions.  The  main  difference  between  the  two  theories 
is  that  the  Onsager  theory  uses  the  second  virial  expan¬ 
sion  of  the  interaction  energy  whereas  the  Flory  theory 
uses  the  lattice  model  to  simulate  the  excluded  volume 
effect.  Because  of  its  virial  expansion,  the  Onsager  theory 
is  more  applicable  to  dilute  solutions.  Introducing  higher 
virial  terms  may  extend  the  Onsager  theory  to  concen¬ 
trated  solutions.59, 60  In  contrast,  the  Flory  theory  makes 
some  assumptions  valid  only  in  a  perfectly  ordered  phase. 
These  assumptions  restrict  the  Flory  theory  to  high  con¬ 
centrations  and  a  highly  ordered  phase.  Both  the  Onsager 
and  Flory  theories  predict  an  isotropic-nematic  phase  tran¬ 
sition  when  the  polymer  concentration  is  sufficiently  high 
(or  equivalently  when  the  temperature  is  sufficiently  low). 

The  Leslie-Ericksen  (LE)  vector  theory24, 73  is  one  of  the 
earliest  constitutive  theories  of  liquid-crystals.  It  extends 
Oseen  and  Frank’s  continuum  theory  of  elasticity  of  liquid 
crystals  at  rest  to  model  flow  behavior.44  The  LE  model 
takes  the  microstructure  of  the  material  into  account  by 
postulating  a  unit  “director”  that  characterizes  the  direction 
of  local  ordering  of  the  polymer  rods.  The  theory  assumes 
a  uniquely  defined  principal  axis  exists  everywhere.  The 


LE  model  has  been  successfully  applied  to  studying  the 
flow  of  nematic  rod  dispersions  of  low  molar-mass.20 
Since  it  assumes  the  linear  dependence  of  the  stress  tensor 
on  the  deformation  rate,  the  LE  theory  is  limited  to  low 
deformation  rates  and  cannot  describe  nonlinear  rheologi¬ 
cal  properties,  such  as  the  shear  thinning  of  the  viscosity 
of  LCPs.  Since  the  theory  also  assumes  an  identifiable 
principal  direction  of  orientation  everywhere,  it  must  be 
extended  to  resolve  singularities  which  must  arise  in  the 
core  of  topological  defects.  This  work  has  a  rich  history 
which  the  reader  can  find  in  some  excellent  manuscripts 
on  the  subject.63,64, 103 

Kinetic  theories  for  liquid-crystalline  polymers  based 
on  the  dynamics  of  the  nonequilibrium  orientational  distri¬ 
bution  function  have  been  developed  by  Siegfried  Hess47 
and  Masao  Doi.23  The  basic  idea  behind  the  Doi-Hess 
theory  is  to  treat  the  LCP  as  a  suspension  of  rigid  rod¬ 
like  nematogenic  molecules  and  describe  the  ensemble 
with  an  orientational  probability  density  function  (pdf). 
The  orientational  pdf  evolves  according  to  a  nonlinear 
Smoluchowski  (Fokker-Planck)  equation  where  the  forces 
acting  on  the  rodlike  molecules  include  hydrodynamic, 
Brownian  and  intermolecular  forces.  The  orientational 
order  present  in  a  nematic  LCP  derives  from  intermolec¬ 
ular  forces  of  various  origins.  Such  molecular  interactions 
are  represented  in  the  Smoluchowski  equation  by  a  mean- 
field  potential  known  as  the  nematic  potential.  The  two 
commonly  used  intermolecular  potentials  are  the  Onsager 
potential  and  the  Maier-Saupe  potential.  The  Onsager 
potential  is  based  on  excluded-volume  interactions  of  a 
single  rigid  rodlike  test  molecule  with  the  other  surround¬ 
ing  molecules.  The  Onsager  potential  can  be  written  as 

VonsagerM  =  VKT  f  /3(m,  m')p(m')  dm'  (1) 

J\m'  1=1 

where  m  and  m'  are  unit  vectors,  p(m)  denotes  the 
probability  density  function  that  a  random  molecule  in 
the  ensemble  has  orientation  m  at  time  t,  /3(m,  m')  = 
2dL2\m  x  m'|  represents  the  excluded-volume  interaction 
between  two  rodlike  molecules  with  orientations  m  and  in', 
v  is  the  number  density  of  the  rigid  rodlike  molecules  with 
diameter  d  and  length  L.  kB  is  the  Boltzmann  constant,  and 
T  is  temperature.  One  can  approximate  the  magnitude  of 
the  cross  product  by  the  following  expression: 

|mxmj  =  V  sin2  0  =  Vl  —  cos2  0 

«  1  —  ^  cos2  0  =  1  —  ^  (m  •  m')2  (2) 

where  ©  is  the  angle  between  m  and  m'.  Using  the  above 
approximation,  one  can  simplify  the  Onsager  potential  to 
the  Maier-Saupe  potential: 

3  UkBT 

Vms  = - 2 — mm  :  M  (3) 
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where  U  is  the  nematic  strength,  and  M  is  the  second 
moment  of  the  orientation  distribution 

M  =  (mm)  =  /  mmp(m)  dm  (4) 

■'|m|=l 

The  Maier-Saupe  potential  (3)  was  originally  developed 
for  studying  phase  transitions  in  thermotropic  liquid  crys¬ 
tals  and  it  has  been  widely  adopted  in  the  study  of  liquid 
crystalline  polymers.  The  Maier-Saupe  is  analytically  less 
complicated  in  that  it  is  completely  specified  by  the  sec¬ 
ond  moment.  In  contrast,  the  Onsager  potential  depends 
on  the  whole  pdf. 

The  fact  that  the  Maier-Saupe  potential  is  completely 
specified  by  the  second  moment  also  inspires  a  further  sim¬ 
plification  in  which  the  Smoluchowski  equation  is  approx¬ 
imated  by  an  equation  for  the  second  moment.  As  a  result, 
the  macroscopic  properties  of  interest  at  any  material  point 
depend  not  on  the  distribution  function  p  directly,  but 
rather  on  low  moments  of  p.  Indeed,  this  projection  is  how 
one  makes  contact  with  the  Landau-deGennes  approach 
which  is  to  posit  closed  equations,  the  free  energy,  and 
stress  constitutive  law  explicitly  in  terms  of  the  second 
moment  of  the  pdf.  Below  we  review  briefly  the  approach 
of  closure  approximation.  This  approach  has  been  widely 
adopted  in  the  development  of  numerical  methods  for  solv¬ 
ing  the  closure -based  class  of  problems  for  dilute  polymer 
solutions  (e.g.,  “Boger  fluids”)  and  other  complex  fluids 
including  liquid  crystal  polymers.  For  simplicity,  here  we 
restrict  ourselves  to  the  two-dimensional  Doi-Hess  theory 
in  the  absence  of  flow. 

We  parametrize  the  rod  axis  by  m  =  {mx,  m2).  Let  9  be 
the  angular  coordinate  of  m.  We  have 

m  j  =  cos  9,  m2  =  sin  9 

dml  dm2  (5) 

- =  —  m2,  - =  m , 

99  2  99  1 

The  stochastic  evolution  of  one  polymer  rod  is  described 
by  a  Langevin  equation:70 

Y'(O') 

d0  =  -D——!-dt  +  V2D-dW  (6) 

kBT 


where  D  is  the  diffusion  constant,  V(9)  is  the  total 
potential  including  the  intermolecular  potential  and  other 
potentials  in  the  absence  of  external  fields,  and  W(t)  is 
the  Wiener  process  (Brownian  motion).  In  general,  the 
Langevin  equation  approach  is  limited  because  of  the  non¬ 
locality  of  the  excluded  volume  potential,  but  more  so 
because  there  are  many  degeneracies  due  to  rotational 
invariance  of  equilibria,  and  in  shear  flow  there  are  many 
limit  cycles  and  bifurcations  which  are,  thus  far,  not 
detectable  from  statistics  of  stochastic  differential  equation 
paths.  One  popular  approach  is  to  use  a  mesoscopic  prob¬ 
ability  density  function  to  describe  the  molecules  and  then 
use  the  Maier-Saupe  potential  which  gives  the  interaction 


potential  between  a  (test)  polymer  rod  and  the  surrounding 
polymer  rods  described  by  mesoscopic  probability  density. 
The  Maier-Saupe  potential  can  be  expressed  as 

3  UknT  (  \\ 

VMS  = - 2 — mm:  ^M—  -J  =  —2NkBTmm:Q  (7) 


where  N  —  3U/4,  and  Q  is  the  orientation  tensor  (traceless 
normalization  of  the  second  moment) 


Q  =  (mm) 


1 

2 


0 11  012 
021  022 


(8) 


The  orientation  tensor  Q  is  symmetric  and  traceless,  so 

012  —  021  and  011  =  —  022- 

The  probability  density  function  p(9,  t)  is  described  by 
the  Smoluchowski  equation 


s_p=  UY^iO)  sp) 

9t  r99\  kBT  r  99 J 


(9) 


where  Dr  is  the  rotational  diffusion  coefficient  of  polymer 
rods.  For  simplicity  we  set  Dr  =  1,  which  is  equivalent 
to  a  re-scaling  of  time  for  constant  Dr ;  the  more  general 
orientation-dependent  diffusion  constant  does  not  lead  to 
qualitatively  different  results.  To  derive  an  equation  for  the 
orientation  tensor,  we  first  use  (5)  to  derive  the  identity: 


—  (mm)  —  (mm):Q 


99  \99 


=  (-4(mm)((mm)  :  Q)  +  4(mm)Q)  (10) 


The  proof  of  this  equation  is  straightforward.  We  first  show 
that  the  (1,  1) -component  is  the  same  for  both  sides  even 
before  taking  the  average. 

=  — 2m1m2[— 2m1m2(<211  -  Q22)  +  2 (m\  -  m2)g12] 

=  — 4w([(— m2)0n  +  m2<222]  —  4m1m2(2m(  —  1)012 
=  -4  m\[m\Qn  +  m22Q22]  +4m]Qn 
-  4 m\ [2m1m2Ql2]  +  4m lm2Ql2 
=  — 4{(nun)[(mm):Q]}11+4{(mm):Q}11  (11) 

Next  we  check  the  (1, 2)-component. 

=  <('«i  -  m\)[-2mlm2{Qn  -  022)  +  2{m\  -  »?;)012]) 
=  (-4mlm2(m\-m22)Qn+2(m\-ml)2Ql2)  (12) 

On  the  other  hand. 


9  , 

9  ,  x  ^ 

(  — (mm) 
\99y 

—  (mm)  :  Q 
99 

9 

9 

—  (mm) 
99 

—  (mm)  :  Q 
99y 

(— 4(mm)[(mm)  :  Q]  +  4(mm) :  Q)1>2 
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=  (— 4(m1m2)(mi2i1  +m\Q21  +  2mlm2Qn) 

+  A{m\Qn  +  mlm2Q22)) 

=  {-Ampn^m]  -  m22)Qn  +  2{m\  -  m22)Ql2) 

+  {2{2m\  -  l)Qn  +  Amlm2Q22) 

=  -  m2)Qu  +  2 (m\  -  m22)Q12 } 

+  ( 2(ml-ml)Ql2  +  Am1m2Q21 ) 

=  (— 4m1m2(Wj  -  m2)Qn  +  2{m]  -  m22)Qn } 

+  {2{m\  —  m2))Ql2  +  (Amlm2)Q22 
=  {—Ampn2(m\  -  m22)Qn  +  2{m]  -  m22)Ql2) 
+2(Gn  —  622)612  +  4<212(322 

=  {—Ampn2{m\  -  m2)6n  +  2(m2  -  m22)Ql2) 

—  4622612  4"  4612622 

=  (— 4m1m2(mj  —  m2)6n  +  2{m\  —  m2)Ql2)  (13) 


This  completes  the  proof. 

Combining  the  identity  (10)  with  Smoluchowski  equa¬ 
tion  (9),  we  find 


dQ 

dt 


/  (21  — 4mm  )pdQ 

J  |m|=l 

+  2 N  [  (— 4(mm)[(mm)  :Q]+4(mm)-Q)pr/0 

^|m|=l 

21  —  4(mm)  —  8Al(mmmm)  :  Q  +  81V(mm)  •  Q 
— 4Q  +  8/vYq+  ^  Q  —  8M(mmmm)  :  Q  (14) 


Therefore,  we  have  derived 


Orientation  tensor  model: 

—  =  -4Q  +  8Al{  Q+-)Q-8Af(mmmm)  :Q  (15) 
dt  \  2) 

In  (15),  the  evolution  of  Q  is  affected  by  the  fourth 
moment  (mmmm).  The  coupling  between  moments  in 
(15)  requires  some  closure  approximation.  The  quadratic 
Doi  closure  rule  assumes 


Doi  closure: 

(mmmm)  :  (•)  =  (mm)  (mm)  :  (•) 

=  MM  :  (•)  =  (M  :  (*))M 


(16) 


Different  closure  rules,  such  as  Hinch-Leal,52  Tsuji-Rey,94 
have  been  studied.  We  list  some  2D  closure  rules  below. 
Tsuji-Rey  closure: 


(mmmm)  :(•)  =  - 


(Q:(*))Q+Q2-(*) 

+  Q.(.)-Q+(.).Q2 
—  ((Q  •(.)):  Q)I+^(Q :(.))! 


(17) 


Hinch-Leal  1  closure: 


(mmmm) :  (•)  =  ^  [6M  •  (•)  •  M  -  (•) :  (MM)  ^ 


-2((.):MM-(.):M)I] 


Hinch-Leal  2  closure: 


(mmmm)  :  (•)  =  M((»)  :  M) 


+  2[M  •  (•)  •  M  -  M2  ((•) :  M2)/(I :  M2)] 


-a(M) 


[(.)-M 

315 v  21 LV  ’ 

+  M  •  (•)  —  ((•) :  M)I] 


[  where  a(M)  =  exp[2(I  —  3M2  : 1)/(I  —  (•)  :  M2)I] 

(19) 

Detailed  studies  on  closure  rules  can  be  found  in  Ref.  [28]. 
As  illustrated  here,  the  orientation  tensor  theory  is  either 
derived  or  recovered  from  the  kinetic  theory  for  the  density 
function  by  projection  onto  a  second-moment  description 
using  closure  rules. 

To  give  a  more  mathematical  description  on  the  ori¬ 
entation  tensor  Q,  let  us  temporarily  leave  2D  rods  in 
favor  of  3D  rods.  The  orientation  tensor  Q  is  the  basis 
for  micron-scale  light  scattering  measurements  of  primary 
axes  (directors),  degrees  of  molecular  alignment  (birefrin¬ 
gence),  and  normal  and  shear  stress  measurements.  The 
spectral  representation  of  Q  gives  mesoscopic  directors 
and  order  parameters.  In  three-dimensional  space  (3D)  the 
eigenvalues  of  Q,  d:  —  1  /3,  are  the  order  parameters,  which 
describe  the  averaged  degrees  of  orientation  with  respect 
to  the  eigenvectors  of  Q,  n(,  which  in  turn  are  the  meso¬ 
scopic  principal  optical  axes.  The  ellipsoid  of  Figure  2 
depicts  a  geometric  representation  of  the  orientation  tensor 
Q  in  3D  space. 

An  experimentalist  can  probe  the  optical  axes  by  vary¬ 
ing  the  plane  transverse  to  the  light  signal;  the  principal 
axes  corresponds  to  a  plane  with  local  maximal  birefrin¬ 
gence  (or  anisotropy).  Birefringence  in  the  plane  of  n,,  n; 
(or  any  plane)  is  measured  by  differences  in  the  principal 
degrees  of  orientation,  Blj  —  \dl  —  d} ].  We  further  note  that 
df  >1/3  implies  that  on  an  average  the  molecules  align 
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Fig.  2.  A  geometric  representation  of  the  mesoscale  3D  orientation  ten¬ 
sor  Q.  The  axes  lengths  d,’s  are  given  by  eigenvalues  (dt  —  1/3)  of  the 
symmetric  matrix  Q,  whereas  the  axes  are  determined  by  the  orthonormal 
frame  of  eigenvectors  of  Q. 

toward  the  principal  axis  n,,  whereas  dt  <1/3  implies 
greater  alignment  with  the  plane  perpendicular  to  n,.  The 
principal  axis  associated  with  the  largest  value  of  d,  is 
usually  referred  as  the  “major  director.”  Mathematically, 
a  spectral  representation  of  Q  is  given  by 


X>«=  1 
1=1 

dj  =  (n,m)  =  (cos2(n,m)) 

where  (n,m)  is  the  angle  between  the  directors  n,  and  m. 
The  isotropic  state  corresponds  to  the  case  where  d  {  = 
d2  =  d3=  1/3  or  equivalently  Q  =  0,  and  the  ellipsoid  of 
Figure  2  reverts  to  a  sphere  with  random  mesoscale  order. 

In  the  presence  of  flow  field,  the  probability  density 
function  p(m,  x,  t)  is  governed  by  the  Smoluchowski 
equation  below,  in  which  the  time  has  been  re-scaled  with 
the  rotational  diffusion  coefficient. 

Smoluchowski  equation: 

=  SR  ■  (pSRpS)  —  SR  ■  (m  x  mp)  (20) 

In  (20),  D/Dt  —  3 /fit  +  v  ■  V  denotes  the  material  time 
derivative,  p,(m)  =  lnp(Tn)+  V(m)/kHl'  is  the  normalized 
chemical  potential  density,  SR  =  m  x  d/d m  is  the  rotational 
gradient  operator,  and  rh  is  the  Jeffery  orbit  of  ellipsoids 
given  by 

m  =  Om  +  fl(Dm-D:mmm)  (21) 

In  (21),  D  and  fl  are  the  rate-of-strain  and  vorticity  ten¬ 
sors,  and  a  is  the  aspect  ratio  parameter  defined  as  a  = 
(r2+l)/(r2— 1)  where  r  is  the  ratio  of  the  rod  length  or 
plate  thickness  to  the  diameter.  0  <  a  <  1  corresponds  to 
the  case  of  rod-like  molecules  and  —  1  <  a  <  0  for  the  case 
of  platelets.  In  the  presence  of  external  electric/magnetic 
fields,  the  potential  V(m)  in  /i(m)  above  will  include 
terms  in  addition  to  the  Maier-Saupe  potential. 

The  development  of  the  Doi-Hess  theory  has  given  a 
large  impetus  to  the  study  of  rheology  of  liquid-crystalline 


polymers.  For  planar  linear  flows  the  monodomain  prob¬ 
lem  has  been  solved.37  What  remains  open  is  full  3D  linear 
flows. 

2.  MATHEMATICAL  STUDIES  OF  PURE 
NEMATICS 

The  isotropic-nematic  (I-N)  phase  transition  in  hard  rod 
gases  and  liquids  is  a  classical  topic,  which  was  first 
explained  theoretically  by  Onsager  in  terms  of  excluded 
volume  potential.87  Onsager’s  approach  was  based  on  a 
variational  model  where  the  second  virial  coefficient  was 
calculated  for  excluded  volume  interaction  between  the 
rods.  His  calculation  resulted  in  a  non-analytic  |m  x  m'| 
potential  where  m  represents  the  direction  of  the  rod.  Later 
Maier  and  Saupe  (3)  re-examined  the  I-N  transition  with 
a  simpler  potential  in  the  form  of  —  (m  •  m')2  that  now 
bears  their  names.80  The  Maier-Saupe  potential  is  a  quar- 
tic  approximation  of  the  Onsager  potential,  which  affords 
sufficient  degrees  of  freedom  to  capture  the  hysteresis 
loop  in  an  order  parameter  representation  of  isotropic  and 
anisotropic  equilibria.  The  interaction  employed  in  the 
Maier-Saupe  theory  of  the  nematic  state  has  been  gener¬ 
alized  in  a  manner  consistent  with  the  asymmetry  of  the 
molecules  that  exhibit  a  nematic  phase.45  The  bifurcation 
in  Onsager’s  model  of  the  isotropic-nematic  transition  in  a 
three-dimensional  system  of  hard  rods  as  well  as  the  corre¬ 
sponding  two-dimensional  system  of  hard  lines  was  inves¬ 
tigated  in  Ref.  [58].  It  was  found  that  the  existence  and 
order  of  a  phase  transition  depend  on  both  the  direction 
of  bifurcation  and  on  properties  of  the  global  solutions. 
In  1986  Doi  and  Edwards23  employed  a  second-moment 
closure  approximation  for  the  Smoluchowski  equation  for 
the  probability  density  function  (PDF)  to  further  illustrate 
robustness  of  the  I-N  transition  with  coarse-grained  mod¬ 
els  and  an  excluded- volume  potential.  In  Ref.  [40]  a  simple 
proof  was  given  to  show  that  all  steady  states  of  the  Doi- 
Hess  tensor  model  are  unaxial.  In  fact,  the  homogeneous 
nematic  dynamics  of  liquid  crystalline  polymers  without 
flow  is  described  by  the  following  Doi-Hess  orientation 
tensor  model: 

jfQ=-F(Q)  (22) 

Here 

F(Q) =  ( 1  -  j  )  Q -  N(Q ■  Q) +  N(Q :  Q)  (q+  0  (23) 

where  Q  is  a  second-order,  symmetric,  traceless  tensor  that 
describes  the  average  molecular  orientation  with  respect  to 
a  probability  distribution,  N  is  the  dimensionless  polymer 
concentration  that  dictates  the  strength  of  the  Maier-Saupe 
potential.  The  orientation  tensor  Q  is  said  to  have  uniaxial 
symmetry  if  two  of  its  eigenvalues  are  identical;  otherwise, 
it  is  said  to  exhibit  biaxial  symmetry.  In  order  to  show  that 
all  homogeneous  equilibria  of  (22)  are  uniaxial,  we  first 
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diagonalize  the  equation  F(Q)  =  0  using  orthogonal  frame 
of  eigenvectors  of  Q.  Then  we  show  that  all  three  eigen¬ 
values  of  Q  satisfy  a  quadratic  equation;  therefore  there 
are  at  most  two  distinct  eigenvalues  out  of  three,  which 
is  the  definition  of  uniaxiality.  However,  the  dynamics  of 
the  orientation  tensor  Q  is  not  necessarily  uniaxial.  The 
stability  of  each  steady  state  was  examined  in  full  system 
in  Ref.  [40]. 

The  proof  that  all  kinetic  model  equilibria  are  uniaxial 
is  much  harder  and  remained  open  until  quite  recently, 
when  mathematical  analysts,  led  by  Constantin,  Kevrekidis 
and  Titi, 13-15  began  to  revisit  Onsager’s  seminal  papers  and 
provide  modern  rigorous  proofs  about  stationary  solutions. 

In  Refs.  [14, 15,26,77]  the  authors  conducted  a  detailed 
study  of  the  2-D  model  where  the  rod  orientation  is  con¬ 
fined  to  a  unit  circle.  In  Ref.  [13]  Constantin  et  al.  also 
considered  the  3-D  model  where  the  rod  orientation  lies 
on  a  unit  sphere.  They  provided  an  upper  bound  on  the 
number  of  critical  points  from  the  analysis  of  a  transcen¬ 
dental  equation  for  the  orientational  order  tensor.  How¬ 
ever,  a  complete  classification  of  all  equilibria  was  not 
addressed.  Incidentally,  it  was  in  2005  that  three  differ¬ 
ent  proofs  appeared  to  show  that  all  anisotropic  equilib¬ 
rium  solutions  of  the  Smoluchowski  equation  must  be 
axial  symmetric  on  a  sphere  with  Maier-Saupe  interaction 
potential.25, 75, 115  In  both  Refs.  [25  and  75]  the  proof  of  the 
axial  symmetry  was  based  on  spherical  harmonics  expan¬ 
sions.  Stability  of  critical  points  with  Maier-Saupe  inter¬ 
action  potential  was  also  considered  in  Ref.  [25]  without 
details.  Very  recently  in  Ref.  [106]  the  axial  symmetry  of 
several  different  cases,  including  pure  nematic  polymers, 
dipolar  nematic  polymers  and  polymers  in  higher  dimen¬ 
sional  space,  was  unified  by  a  key  inequality.  In  the  fol¬ 
lowing  we  are  going  to  briefly  summarize  these  different 
proofs.  For  more  details,  see  the  Refs.  [25,75, 106, 115]. 

In  Ref.  [25]  Fatkullin  and  Slastikov  provided  a  rig¬ 
orous  analysis  and  classification  of  all  critical  points  of 
Onsager’s  functional.  The  Onsager  free  energy  functional 
can  be  expressed  as  the  sum  of  an  entropic  term  and  an 
interaction  term  and  it  has  the  expression 


Tp(s)  lnp(i) 


a  saddle-node  bifurcation  occurs.  The  important  tool  used 
in  Ref.  [25]  is  the  spherical  harmonics.  Similar  approach 
has  been  used  in  Ref.  [75]. 

In  Ref.  [115]  we  provided  a  proof  based  on  elemen¬ 
tary  calculus.  In  the  absence  of  external  field  and  flow, 
the  governing  equation  for  stiff  rod  nematic  polymers 
interacting  by  the  Maier-Saupe  potential  is  given  by  the 
Smoluchowski  (or  Fokker-Planck)  equation: 


dt  du  \  kB  T  d ii  du  J 

V(u,  [p])  =  —  bkBT(u® u)  :  (u® u)  (25) 


(u<g>u)  =  /  (u<g>u)p(u)t/u 
Js 

where  p  is  the  orientational  distribution  function,  u  a  unit 
vector  specifying  the  orientation  of  a  rodlike  molecule, 
D  the  rotational  diffusion  coefficient,  V  the  Maier-Saupe 
short-range  nematic  potential,  (u  <g>  u)  the  second  moment 
of  p,  kB  the  Boltzmann  constant,  T  the  absolute  tempera¬ 
ture,  and  b  the  strength  of  the  nematic  potential.  In  equi¬ 
librium  the  solution  of  the  Smoluchowski  equation  is  given 
by  the  Boltzmann  distribution: 


p(u)  = 


exp(— V(u,  [p])) 

Is  exp(— V (u,  [p]))  du 


(26) 


Since  (u  <8>  u)  is  symmetric,  we  choose  a  Cartesian  system 
such  that  (u<S>u)  is  diagonal: 


(u®u)  = 


0 


V 


o 


0 


0  \ 


0 


0 


(27) 


where  u  =  (tq,  u2,  iq)  and  r;  =  (w2)  —  1/3.  If  rl  —  r2  — 
r3  =  0,  then  the  system  is  in  isotropic  state;  If  there  are  two 
distinct  eigenvalues  among  r, ,  r2  and  r3 ,  then  the  system 
is  in  axisymmetric  (uniaxial)  state.  Upon  diagonalization 
of  the  second  moment  (u<S>u),  the  probability  distribution 
function  takes  a  simple  form: 


p(s)  J  U(s,  s')p(s')  ds' 


ds  (24) 


P(u) 


exp|  k(r,  u2  +  r2u\  +  r2u\  +  1/3)] 
fs  exp[k(r,  u2  +  r2u\  +  r2u\  +  1/3)]  du 


(28) 


where  p(s)  is  a  probability  density  function  describing 
rod  orientations  for  a  system  of  liquid-crystalline  suspen¬ 
sions,  s  is  the  orientation  parameter  which  lies  on  a  unit 
sphere  S2,  T  is  the  temperature  (and  the  Boltzmann  con¬ 
stant  kB  is  set  to  unity),  and  U ( s ,  s')  is  the  Maier-Saupe 
potential.  The  uniform  density  p  =  1/477  describes  the 
isotropic  phase  where  the  molecules  are  randomly  ori¬ 
ented.  The  main  result  in  Ref.  [25]  is  that  the  uniform 
density  is  a  unique  critical  point  when  T  >  T*  sa  0.149 
and  an  isolated  minimizer  for  T  >  Tc  =  2/15.  At  T  =  T* 


where  b  is  treated  as  a  variable  for  mathematical  conve¬ 
nience.  Then  we  introduce  three  random  variables: 

hi  =  m2-(«2),  i  =1,2, 3  (29) 

and  show  that  they  have  following  properties: 

1.  <A,->  =  0,  i  =  1,2, 3;  EL  */  =  <). 

2.  Suppose  ;■]  >  r2.  Then  for  b  >  0,  we  have 
(k3(/q  —  h2))  <  0.  Note  that  (/z3/q)  is  the  correlation  of 
u\  and  u\  while  (h3h2)  is  the  correlation  of  u 2  and  u2. 
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We  further  define  functions 

Ri(b)  =  i=  1,2,3  (30) 

which  satisfy 

1.  R,(0)  =  0,  i  =  1,2,3 

2.  If  (r, ,  r2,  r3)  corresponds  to  an  equilibrium  for  b  =  b{), 
then  Rj(b0)  —  rh  i  =  1,  2,  3 

We  prove  the  axisymmetry  of  the  equilibrium  solutions  by 
contradiction.  Suppose  there  is  a  set  of  distinct  r, ,  r2,  r3  for 
b  —  b0>0.  Without  loss  of  generality,  assume  r,  >  r2>  r3. 
Because  r,  +  r2  +  r3  —  0,  we  have  rx  >  0  and  r3  <  0.  Con¬ 
sider  a  function 


F{b)  =  r2Rl{b)-rlR2{b)  (31) 

It  is  clear  that  F(0)  =  F(b0)  —  0.  To  establish  contradic¬ 
tion,  we  carried  out  detailed  analysis  in  Ref.  [115]  to  show 
that 


2F\b)  =  (c  —  r2)(r2  —  r3)(h2(h3  —  h{)) 

+  {r2  -  r3)(-r3)(h3(h2  -  hx)) 

+  (.rl-r2)(rl)(h2(h3-h2))  >0 

for  b  >  0  (32) 

since  all  the  coefficients  and  correlation  terms  are  positive. 
(32)  implies  that  F{b )  is  a  monotonically  increasing  func¬ 
tion,  which  contradicts  with  the  fact  F’(O)  =  F(b0)  =  0. 

To  estabilish  the  number  of  equilibria,  we  select  the  axis 
of  symmetry  as  the  z-axis  and  use  the  spherical  coordinates 
to  denote  the  unit  vector  u: 

u  =  (ul,u2,  m3)  =  (sin  <p cos  9,  sin  </>  sin  9,  cos  <f>)  (33) 

Then  the  Maier-Saupe  potential  has  the  expression 

V(0)  =  —  r  cos2  cf>  (34) 

where  r  =  {b/ 2)(3(cos2  <fi)  —  1)  is  the  Flory  order  param¬ 
eter  and  satisfies  the  equation 

r[l-bf(r)]=0  (35) 

Here  the  function  f(r)  is  given  by 

=  /o*  fe>2(l  ~  tu2)  exp(rw2)  da> 

f0‘  exp (rco2)  dco 

and  can  be  shown  to  have  the  following  properties: 

I-  ./  (0)  —  2/15 

2.  0<  /(/-)<  1/4 

3.  liml.^_00/(r)  =  0,  \\mr^x  f(r)  =  0 

4.  There  exists  r*  >  0,  such  that  =  0,  f'(r)  >  0  for 
r  <  r*  and  /'(r)  <  0  for  r  >  r* 

Then  the  conclusion  on  the  number  of  equilibria  follows 
immediately.  More  specifically, 


1.  For  b  <  b*,  there  exists  one  solution  r  —  0. 

2.  For  b  =  /;*,  there  are  two  solutions:  r  =  0  and  r  =  r*  — 
2.17828797  •  •  • . 

3.  For  b*  <  b  <  15/2,  there  are  three  solutions  0,  r,  and 
r2  where  i\  is  between  0  and  r*,  and  r2  is  bigger  than  r* . 

4.  For  b  =  15/2,  there  are  two  solutions  0  and  r2  where 
r2  is  bigger  than  r*. 

5.  For  b  >  15/2,  there  are  three  solutions  0,  r,  and  r2 
where  rx  is  negative  and  r2  is  bigger  than  r*. 

Here  b*  =  \/f{r*)  =  6.73 1486  •  •  • . 

In  contrast,  in  the  Doi-Hess  orientation  tensor  theory 
the  uniaxial  equilibrium  branches  as  functions  of  the  con¬ 
centration  parameter  N  (inversely  proportional  to  the  tem¬ 
perature)  are  given  by:40 


■so  =  0. 


l±3Vl-8/3  N 


(37) 


Figure  3  depicts  the  isotropic-nematic  phase  diagram 
for  the  Maier-Saupe  interactions  with  the  Doi-Hess  ori¬ 
entation  tensor  theory.  It  was  further  shown  in  Ref.  [40] 
that  the  “oblate  phase”  (where  the  order  parameter  s  <  0) 
is  stable  to  uniaxial  perturbations,  but  unstable  to  biax¬ 
ial  perturbations.  Thus  it  is  unstable.  This  explains  why 
“oblate  phase”  is  not  observed  in  equilibrium  experiment; 
the  oblate  phase  turns  out  to  play  a  prominent  role  in  shear 
flows  when  there  are  different  limit  cycles  (tumbling  and 
wagging,  for  example)  at  different  heights  in  the  shear  gap. 

Rigorous  analysis  for  steady  state  solutions  with  the 
Onsager  potential  is  still  unavailable.  The  main  difficulty 
lies  in  the  fact  that  the  Onsager  potential  depends  on 
the  whole  function  of  the  probability  density  whereas  the 
Maier-Saupe  potential  depends  only  on  the  second  moment. 


Fig.  3.  Isotropic-nematic  phase  diagram  for  the  Maier-Saupe  inter- 
molecular  potential  with  the  Doi  closure  tensor  theory.  The  horizontal 
axis  represents  the  strength  of  the  nematic  potential  whereas  the  vertical 
axis  denotes  the  nematic  order  parameter. 
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For  extended  nematics  which  consists  of  an  ensemble  of 
rigid  rods  with  inherent  dipoles,  the  Smoluchowski  equa¬ 
tion  can  be  written  as 

dt  <9u  \ksT  duH  <9u /  (38) 

U( u,  [p])  =  —akBT( u)  u  —  M'Br(u<8>u)  :  (u®u) 

Here  a  is  strength  of  the  dipole-dipole  interaction.  In 
Ref.  [116]  it  was  proved  that  all  stable  equilibria  of 
rigid,  dipolar  rod  dispersions  are  either  isotropic  or  prolate 
uniaxial.  It  was  further  demonstrated  that  unstable  non- 
axisymmetric  equilibria  exist  and  therefore  the  stability  is 
essential  in  establishing  the  axisymmetry.  A  recent  work 
by  Lee  et  al.72  studied  dipole-induced  first-order  phase 
transitions  of  nano-rod  monolayers  where  the  total  poten¬ 
tial  for  dipolar  Brownian  rods  in  equilibrium  is  given  by 
(38).  It  was  shown  that  the  coupling  of  a  dipole  poten¬ 
tial  to  excluded  volume  is  sufficient  to  re-instate  a  first- 
order  phase  transition  of  rods  confined  to  two-dimensions. 
Extending  the  work  to  imposed  elongational  or  shear  flow, 
or  a  magnetic  field  remains  to  be  explored. 

3.  EFFECTS  OF  FLOW  FIELDS 

Even  though  nematic  liquid  crystalline  polymers  (LCPs) 
have  been  studied  for  quite  a  long  time,  their  technologi¬ 
cal  applications  remain  dominated  by  high  tensile  strength 
fibers,  such  as  Kevlar.  This  is  mainly  due  to  the  intrinsic 
difficulty  of  maintaining  a  monodomain  orientational  state 
at  the  macroscopic  level.  Understanding  the  microstructural 
orientation  state  in  any  processing  flow  other  than  fiber 
spinning  becomes  extremely  important  in  any  future  devel¬ 
opments  of  LCPs  as  structural  materials.  Many  numeri¬ 
cal  studies  have  been  performed  on  the  characterization 
of  the  shear  flow  behavior  of  LCPs  based  on  various 
modeles,  including  the  Leslie-Erickson  theory,50,51  the  Doi 
model,27-29, 66’ 98  and  the  Smoluchowski  equation.32,33,35,36 
These  studies  aimed  to  reproduce,  explain  and  predict 
experimental  discoveries  of  steady  and  transient  modes. 
“By  a  combination  of  theory  and  experiment,  many  steady 
and  transient  shear-induced,  monodomain  modes  have  been 
catalogued  and  named  primarily  on  the  basis  of  director 
response:  steady  alignment  with  primary  director  either 
in  the  shear  plane  (flow  aligning )  or  along  the  vorticity 
axis  ( logrolling );  in-plane  transient  oscillatory  ( wagging )  or 
rotating  ( tumbling )  director  modes;  and  out-of-plane  tran¬ 
sient  director  modes  (kayaking).  Complicated  dynamics  is 
also  possible.”37  The  effect  of  extensional  flow  and  general 
linear  planar  flows  has  also  been  investigated.31,34,37  The 
proof  of  limit  cycles,  except  in  a  simple  2D  tensor  model, 
remains  open  to  our  knowledge  in  3D  tensor  models  and 
in  either  2D  or  3D  kinetic  equations. 

We  are  now  going  to  impose  physical  boundary  condi¬ 
tions  on  the  oppositely  moving,  steady  phates.  This  will 
necessarily  lead  to  spatial  structure  in  the  orientation  field. 


and  thereby  generate  nonlinear  shear  flow  and  potentially 
higher  dimensional  structures  and  flows.  For  the  Doi-Hess 
mesoscopic  orientation  tensor  model,  the  asymptotic,  one¬ 
dimensional  gap  structures,  along  with  the  flow-gradient 
direction,  in  “slow”  Couette  cells  were  derived  in  Ref.  [42]. 
A  simple  geometric  setup  of  this  study  is  sketched  in 
Figure  4. 

The  scaling  properties  of  steady  flow  molecular  struc¬ 
tures  in  slow  Couette  flows  with  equal  elasticity  constants 
in  Ref.  [42]  were  later  generalized  in  several  ways  to  con¬ 
trast  isotropic  and  anisotropic  elasticity  and  to  compare 
Couette  versus  Poiseuille  flow  and  to  consider  dynamics 
and  stability  of  these  steady  states  within  the  asymptotic 
model  equations.16  In  Ref.  [12]  the  classical  problem  of 
the  viscoelastic  response  of  nematic  liquid  crystal  poly¬ 
mers  to  small  amplitude  oscillatory  shear  was  revisited  by 
applying  a  multiple  time  scale  pertrubation  analysis  to  the 
Doi-Hess  rigid  rod  model.  The  analysis  in  Ref.  [12]  was 
restricted  to  in-plane  orientational  configuration.  General¬ 
ization  of  the  analysis  to  full  tensor  degrees  of  freedom 
leads  to  a  five-dimensional  dynamical  system,  which  is  not 
yet  analytically  solved.  In  Ref.  [71]  monolayer  films  of  liq¬ 
uid  crystalline  polymers  were  modelled  with  a  mesoscopic 
two-dimensional  analogue  of  the  Doi-Hess  tensor  model. 
The  weak- shear  steady  and  unsteady  selection  criteria  for 
2D  nematic  polymers  using  various  second-moment  clo¬ 
sures  was  derived.  A  simple  proof  was  given  based  on 
the  Poincare-Bendixon  Theorem  to  show  that  limit  cycles 
(“tumbling  orbits”)  exist  beyond  the  parameter  boundary 
for  the  steady-unsteady  transition.  Finally,  it  was  revealed 
that  the  shear-perturbed  2D  phase  diagram  is  significantly 
robust  to  closure  approximations  than  the  3D  system. 

The  effect  of  homogeneous  elongational  flows  v  = 
v(—x/2,  —y/2,  z)  (v  >  0:  axially  stretching;  v  <  0:  planar 
stretching)  was  studied  in  Ref.  [41]  using  the  Doi-Hess 
orientation  tensor  model.  All  flow-induced  steady  states 
(both  uniaxial  and  biaxial)  and  their  stability  were  investi¬ 
gated.  The  main  conclusions  were  that  in  axial  stretching, 
biaxial  steady  state  solutions  exist  but  they  are  all  unsta¬ 
ble  and  the  only  stable  solutions  are  uniaxial;  in  planar 


Fig.  4.  A  schematic  diagram  of  the  Couette  cell,  with  no  slip  velocity 
conditions  and  tangential  anchoring  at  the  plates. 
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stretching,  above  a  moderate  concentration,  the  only  stable 
solutions  are  biaxial. 

For  the  Smoluchowski  equation,  theoretical  analysis 
becomes  more  challenging.  The  first  elegant  analysis  of  the 
isotropic-to-nematic  phase  transition  for  two-dimensional 
liquid  crystal  polymers  at  high  shear  rates  was  given  by 
Marrucci  and  Maffettone. 82,83  In  Ref.  [120]  the  dynam¬ 
ics  and  steady  states  of  2-D  nematic  polymers  driven  by 
the  weak  shear  were  studied.  We  briefly  review  the  results 
in  120.  For  polymer  rods  in  the  two-dimensional  space, 
the  orientation  of  a  polymer  rod  can  be  described  by  an 
angle  9  and  the  orientation  direction  can  be  represented 
by  the  unit  vector  u  =  (cos  9,  sin  9).  In  the  presence  of  a 
weak  shear  flow  v  =  y(0,  —  x),  the  orientation  distribution 
function  evolves  according  to  the  Smoluchowski  equation 


dp  d 

—  =  D  — 
dt  dd 


-^+lvsH(0)  +  Ks(0))p+d£ 


VMS(d )  =  —  C/(cos 2(0  —  a))  cos 2(0  —  a) 
Xs//(0)=  \  sin  20 


(39) 


where  e  is  the  ratio  of  the  shear  rate  (y)  and  rotational 
diffusivity  ( D )  and  is  called  the  Peclet  number ,  U  is  the 
normalized  polymer  concentration. 

In  the  2-D  case,  it  is  well-known  that  in  the  absence 
of  flow  the  isotropic-nematic  phase  transition  occurs  at 
U  —  2.  In  the  presence  of  an  imposed  weak  shear  there 
is  a  threshold  ( U0  ^  2.41144646)  for  U:  When  U  <  U0 , 
steady  state  solution  exists;  otherwise  there  is  no  steady 
state.  The  slow  time  evolution  driven  by  the  weak  shear 
can  be  captured  by  the  multi-scale  analysis.  It  was  revealed 
in  Ref.  [120]  that  the  order  parameter  of  the  orientational 
distribution  of  rodlike  molecules  is  time  invariant  while 
the  director  changes  with  time.  More  specifically, 

•  The  director  of  the  orientational  distribution  converges 
to  a  stable  steady  state  position  when  the  strength  of  the 
nematic  potential  U  is  less  than  U0. 

•  The  angular  velocity  of  the  director  is  position- 
dependent  and  is  always  positive  when  the  orientational 
distribution  is  temporarily  periodic  (“tumbling”)  and  can 
never  reach  a  steady  state  when  U  is  larger  than  the 
threshold  UQ. 


Recently  a  comprehensive  study  of  the  2-D  nematic  liq¬ 
uid  crystal  polymers  under  an  arbitrary  shear  was  provided 
in  Ref.  [107].  The  study  shows  that  when  the  concentra¬ 
tion  parameter  U  is  smaller  than  a  critical  value  U0,  the 
stable  branch  and  an  unstable  branch  form  a  fold  and  are 
connected  at  U0  and  it  resolves  the  details  of  the  transi¬ 
tion  occurs  around  the  folding  point  of  the  phase  curve. 
As  the  Peclet  number  increases,  the  stable  branch  and  the 
unstable  branch  of  the  fold  are  peeled  off  from  each  other 
and  are  separated.  When  the  Peclet  number  is  larger  than 
0.746,  the  fold  disappears  completely  and  there  is  only  one 
steady  state  for  each  value  of  U. 


The  effect  of  elongational  perturbations  on  nematic  liq¬ 
uid  crystal  polymers  under  a  weak  shear  based  on  the 
two-dimensional  Smoluchowski  equation  was  investigated 
in  Ref.  [118].  Let  u=  (cos  9,  sin  6)  represent  the  orien¬ 
tation  of  rodlike  molecules,  v  =  y(0,  —x)  be  the  shear 
flow,  xE  =  (y/2 )q(—y,  —x)  be  the  elongational  perturba¬ 
tion.  Then  the  Smoluchowski  equation  for  the  orientational 
distribution  function  can  be  written  as 

(y'rSw+vMSmp+d-f0 

VMS(0)  =  —  [/(cos  2(9  —  a))  cos  2(9  —  a)  (40) 
Vps(0)  =  -|0+^(1  +  ?)  sin2  9 

where  £  is  the  Peclet  number.  Applying  the  multi-scale 
analysis,  we  found  that 

•  When  the  elongational  perturbation  is  very  small,  the 
orientational  probability  density  function  tumbles  period¬ 
ically  only  in  an  intermediate  range  of  polymer  concen¬ 
trations;  outside  this  immediate  range  the  orientational 
probability  density  function  converges  to  a  steady  state  and 
the  tumbling  is  arrested. 

•  When  the  elongational  perturbation  is  increased  further 
(by  about  20%  of  the  shear  rate),  the  intermediate  range 
of  tumbling  disappears  and  the  orientational  probability 
density  function  converges  to  a  steady  state  regardless  of 
the  polymer  concentration. 

These  findings  are  consistent  with  various  earlier  results 
based  on  the  Leslie-Ericksen  theory10  or  analogous  3D 
numerical  simulations.43  It  is  worthwhile  to  point  out  that 
even  though  the  rigor  of  multi-scale  asymptotic  analysis 
has  not  yet  been  established,  the  method  itself  has  become 
a  powerful  mathematical  tool  in  recent  applications.  For 
example,  Choate  and  Forest12  applied  it  to  study  the  vis¬ 
coelastic  response  of  nematic  polymers  to  small  amplitude 
oscillatory  shear;  Vicente  Alonso  et  al.102  used  it  to  inves¬ 
tigate  the  nonlinear  dynamics  of  a  nematic  liquid  crystal 
in  the  presence  of  a  plane  Couette  flow  using  a  Landau- 
de  Gennes  model;  Chillingworth  et  al.11  employed  it  to 
describe  the  geometry  and  dynamics  of  nematic  liquid 
crystal  in  a  uniform  shear  flow. 

In  summary,  the  analytical  studies  of  the  Smoluchowski 
equation  for  nematic  liquid  crystal  polymers  under  shear, 
elongational,  or  other  simple  planar  lineal'  flows  have  been 
limited  to  relatively  simple  two-dimensional  case.  Very 
recently  Zhang  and  Zhang111  provided  a  theoretical  study 
on  stable  equilibria  of  liquid  crystals  without  flow  and 
the  stable  dynamic  states  for  nematic  liquid  crystals  under 
weak  shear  flow  for  the  Smoluchowski  equation.  Semi- 
analytical  results  on  the  Smoluchowski  equation  with  the 
Onsager  potential  are  provided  in  Ref.  [105].  Extensions 
to  general  shear  flow  and  general  linear  planar  flow  remain 
unexplored. 


dp  d 

—  =  D  — 
dt  dd 


10 


J.  Com  put.  Theor.  Nanosci.  7,  1-16,  2010 


Zhou  et  al. 


Mathematical  Studies  and  Simulations  of  Nematic  Liquid  Crystal  Polymers  and  Nanocomposites 


4.  NANOCOMPOSITES  AND  EFFECTIVE 
CONDUCTIVITY 


One  particular  area  of  nanotechnology  research  has  been 
carbon  nanotubes  (CNTs).  Ever  since  their  discovery  in 
1991,  nanotubes  have  emerged  as  stars  of  the  chemistry 
world.  They  are  stronger  than  steel,  lightweight,  and  able 
to  withstand  repeated  bending,  buckling,  and  twisting;  they 
can  conduct  electricity  as  well  as  copper  or  semiconduct 
like  silicon;  and  they  transport  heat  better  than  any  other 
known  material.  CNTs  can  be  thought  of  as  a  sheet  of 
graphite,  which  has  its  two  ends  wrapped  together  to  form 
a  hollow  cylinder  and  has  a  cap  enclosing  the  tube  at  one 
of  the  ends.  CNTs  exhibit  a  high  aspect  ratio  (length-to- 
diameter  ratio),  with  the  length  generally  several  p, m  and 
the  diameter  approximately  0.4  to  100  nm.  Nano-elements 
like  CNTs  are  combined  at  low  (0.1-5%)  volume  fractions 
with  traditional  polymeric  materials  to  enhance  a  diverse 
set  of  targeted  properties.  Notable  examples  include  high 
electrical  conductivity  and  strength  of  CNTs  and  the  bar¬ 
rier  properties  of  nano-clay  platelets.  Nanoparticle  suspen¬ 
sions,  ranging  from  DNA,  viruses  and  carbon  nanotubes 
to  mineral  suspensions  with  rod-or  disk  like  particles  also 
form  liquid  crystal  phases,  as  do  micron-sized  colloidal 
particles  in  solvents. 

Recently  Zheng  et  al.112  applied  homogenization  theory 
to  predict  the  effective  conductivity  tensor,  cre ,  for  uniform 
suspensions  of  ellipsoidal  nano-inclusions  of  rigid  rods 
with  conductivity  cr,  and  volume  fraction  02  in  a  matrix  of 
conductivity  oq.  For  simplicity,  all  ellipsoids  are  assumed 
to  have  the  same  geometry  with  three  semi-axes  a ,  b  and  c 
(a  >  b  —  c  for  spheroidal  rods  and  a  —  b>c  for  spheroidal 
platelets).  The  ellipsoids  orient  due  to  excluded-volume 
interactions  and  flow  according  to  a  orientation  probability 
distribution  that  is  the  central  oject  of  the  Doi-Hess  kinetic 
theory.  From  homogenization  theory  of  composites  with 
spheroidal  inclusions  at  low  volume  fractions,  conductivity 
obeys  elliptic  equations  which  are  virtually  identical  for 
thermal,  electric,  and  dielectric  properties.  The  effective 
conductivity  tensor  in  close  form  is  derived  in  Ref.  [112]: 


cr1'  =  oql-l- crl02(cr2-  oq) 


+ 


(T2  +  arl-((72-crl)La 

_ (tr2-o-i)(l  —3  Lfl) _ 

[(cr2  +  o^)  -  (cr2  -  <tx)L  J[o-j  +  ((cr2  -  cr1)La] 

-O(02) 


M(p) 


(41) 


where  I  is  the  3  by  3  identity  matrix,  La  is  the 
spheroidal  depolarization  factor  depending  on  the  aspect 
ratio  r  (r  —  a/b^>  1,  rods;  r  =  c/a«l,  platelets)  of  the 
molecular  spheroids  through  the  relation 


L„ 


1  /1  +  e  , 

—  In  - 1 

2e  V  1  -  e 


(42) 


Here  p  is  the  orientation  probability  distribution  func¬ 
tion  of  quiescent  or  flowing  anisotropic  macromolecules, 


governed  by  the  Smoluchowski  equation  of  the  Doi-Hess 
kinetic  theory  for  quiescent  or  flowing  nematic  polymers; 
M(p)  is  the  second-moment  of  the  orientation  probability 
density  p.  From  the  closed  form,  the  scaling  properties  of 
enhanced  conductivity  versus  volume  fraction  and  weak 
shear  rate  become  explicit.  Not  surprisingly,  the  effective 
conductivity  tensor  inherits  all  the  features  of  the  orien¬ 
tation  tensor  from  the  isotropic-nematic  phase  transition, 
including  hysteresis,  bi-stability,  and  discontinuous  jumps. 

In  Ref.  [113]  we  combined  two  classical  mathematical 
asymptotic  analysis,  slender  longwave  hydro-dynamics30, 41 
and  homogenization  theory  for  composites,112  to  predict 
elongational  flow  enhancements  in  conductivity  of  LCP 
nano-composites.  We  found  that  elongational  flow  dom¬ 
inates  free  surface  and  thermal  effects  on  electrical  and 
thermal  conductivity  enhancement  and  there  is  no  sacri¬ 
fice  in  these  properties  by  producing  much  larger  radius 
fibers. 

Even  though  the  above  studies  are  restricted  to  bulk 
homogeneous  mesophases  of  nematic  nanocomposites  at 
rest  and  in  weak  shear  flows112  and  in  elongational  flow,113 
they  serve  as  benchmarks  and  lay  the  groundwork  for  future 
extensions  to  heterogeneity,  more  general  flow  rates  and 
flow  type,  and  the  incorporation  of  effects  of  the  “inter¬ 
phase”  between  the  matrix  and  nano-inclusions. 

5.  NUMERICAL  SIMULATIONS 

As  we  have  seen  already,  the  governing  equations  for  liq¬ 
uid  crystal  polymers  (and  other  complex  fluids)  are  non¬ 
linear  and  the  only  feasible  way  of  comparing  model 
predictions  with  experimental  observations  is  to  carry  out 
numerical  simulations.  In  1990  Larson68  carried  out  one 
of  the  earliest  numerical  simulations  based  on  the  kinetic 
Doi-Hess  model  with  the  Onsager  potential.  He  solved 
the  Doi  equation  for  the  three-dimensional  time-dependent 
orientation  distribution  function  for  rodlike  nematic  poly¬ 
mers  by  an  expansion  in  spherical  harmonic  functions.  His 
results  confirmed  two-dimensional  calculations  with  the 
Maier-Saupe  potential82  and  also  predicted  negative  first 
normal  stress  difference,  which  is  in  qualitative  agreement 
with  the  experimentally  measured  values.  Later  in  1991 
Larson  and  Ottinger70  used  two  numerical  techniques  to 
solve  the  Doi  equation  with  the  Onsager  potential  again. 
One  of  their  numerical  techniques  is  based  on  spherical 
harmonic  function  expansions,  and  the  other  is  a  stochas¬ 
tic  simulation  technique  which  basically  rewrites  the  Doi 
diffusion  equation  as  a  stochastic  differential  equation  for 
the  time  evolution  of  the  Markovian  process  and  then  inte¬ 
grates  the  equations  of  motion  for  a  large  ensemble  of 
molecules.  They  predicted  flow-aligning,  log-rolling,  wag¬ 
ging,  tumbling  and  kayaking.  As  noted  earlier,  the  abil¬ 
ity  to  detect  bifurcations  among  these  various  different 
orbits  is  beyond  current  methods  of  stochastic  differential 
equations. 
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In  recent  years  the  spherical  harmonics  expansion 
method  has  gained  popularity  in  solving  the  kinetic  Doi- 
Hess  equation.  For  example,  Grosso,  Keunings,  Crescitelli 
and  Maffettone  used  it  in  their  prediction  of  chaotic 
dynamics  in  sheared  liquid  crystalline  polymers.46  Forest’s 
group31-37  carried  out  more  detailed  monodomain  studies 
using  the  spherical  harmonic  expansions,  both  for  mon¬ 
odomains  and  for  spatially  heterogeneous  nematic  polymer 
flows. 

Another  popular  approach  to  model  liquid  crystalline 
polymers  is  to  solve  the  Doi  tensor  model  instead  of  the 
Smoluchowski  equation.  As  we  mentioned  in  the  Introduc¬ 
tion,  the  Doi  tensor  model  is  an  approximation  of  the  Doi- 
Hess  kinetic  model  with  the  the  Maier-Saupe  potential 
and  the  Doi-Marrucci-Greco  potential.  The  Doi-Marrucci- 
Greco  potential  couples  the  Maier-Saupe  short  range  inter- 
molecular  potential  to  a  long  range  potential  to  account 
for  non-local  elastic  interactions.  The  Doi-Marrucci-Greco 
potential  is  given  by 

®dmg  =  -  \  UkB T  Q +  l—  Aq")  :  mm  (43) 

where  /  is  the  persistence  length  of  the  non-local  dis- 
tortional  elasticity  interaction  and  U  is  a  constant  repre¬ 
senting  the  nematic  strength.  Calculations  based  on  the 
Doi  tensor  model  have  been  extensively  carried  out.  For 
example,  in  Ref.  [66]  the  orientation  tensor  is  restricted 
to  be  a  two-by-two  matrix  and  the  coupled  flow  calcu¬ 
lations  are  one-dimesional,  in  Ref.  [61]  the  behavior  of 
the  Doi-Marrucci-Greco  (DMG)  model  for  nematic  liq¬ 
uid  crystalline  polymers  in  planar  shear  flow  is  investi¬ 
gated  where  the  2D  assumption  of  no  spatial  gradients 
in  the  flow  direction  is  retained  in  the  computations.  It 
was  found  in  Ref.  [61]  that  the  DMG  model  was  able  to 
capture  dynamic  behaviors  in  accordance  with  experimen¬ 
tal  observations69  for  the  Ericksen  number  and  Deborah 
number  cascades.  For  increasing  shear  rates  within  the 
Ericksen  number  cascade,  where  gradient  elasticity  is  the 
dominant  stress  contribution,  the  DMG  model  exhibits  four 
distinct  flow  regimes:  stable  simple  linear  shear  flow  at 
low  shear  rates;  stable  roll  cells  at  intermediate  shear  rates; 
irregular  structure  followed  by  large-stain  disclination  for¬ 
mation  and  irregular  structure  preceded  by  disclination 
formation  at  high  shear  rates.  For  increasing  shear  rates 
within  the  Deborah  number  cascade,  where  the  dominant 
contribution  to  the  stress  is  viscoelasticity,  the  DMG  model 
displays  a  stream-wide  banded  texture,  in  the  absence  of 
roll  cells  and  disclinations,  followed  by  a  shear  alignment 
where  the  mean  orientation  lies  within  the  shear  plane 
throughout  the  monodomain.  Very  recently  Klein  et  al. 
extended  their  calculations61  to  three-dimensional  shear- 
driven  dynamics  to  study  polydomain  textures  and  discli¬ 
nation  loops  in  liquid  crystalline  polymers.62 

Finally,  we  want  to  point  out  that  numerical  simulations 
based  on  the  Leslie-Ericksen  theory  on  liquid  crystals  have 


also  been  carried  out  extensively.  For  example,  Rey  and 
his  co-workers50-51  have  used  this  theory  to  study  bifur¬ 
cation  and  banded  textures  of  nematic  polymers  in  simple 
shear  flows.  Very  recently  they  have  conducted  a  nanoscale 
analysis  with  a  continuum  theory  to  predict  a  new  defect- 
forming  mechanism  for  liquid  crystal  interfaces.108 

6.  CONTROLLABILITY  AND  OBSERVABILITY 
OF  NEMATIC  LIQUID  CRYSTALS 

The  controllability  of  viscoelastic  fields  is  a  fundamental 
concept  that  defines  some  essential  capabilities  and  limi¬ 
tations  of  the  resulting  materials.  In  practice,  the  shape  of 
an  extrudate  can  be  controlled  by  varying  the  size  or  shape 
of  an  orifice,  whereas  the  advance  of  the  free  surface  can 
be  controlled  by  varying  the  inflow  into  a  mold.  There  is 
still  a  big  gap  between  theoretical  work  and  practical  prob¬ 
lems.  An  important  issue  of  controllability  is  the  question 
whether  it  is  possible  to  steer  a  system  to  a  desirable  state 
with  a  given  set  of  control  inputs. 

In  Ref.  [89]  the  controllability  of  flows  of  linear 
viscoelastic  fluids  was  investigated  by  Renardy  using 
multi-mode  Maxwell  models.  The  state  of  the  system  is 
characterized  by  the  velocities  and  the  viscoelastic  stresses; 
the  control  input  is  in  the  form  of  the  body  force.  It  was 
found  that  the  system  is  uncontrollable,  unless  the  ini¬ 
tial  conditions  for  the  stresses  satisfy  a  set  of  constraints. 
For  the  degenerate  case  of  creeping  flow  where  the  den¬ 
sity  is  zero  in  the  equation  of  motion,  there  is  no  con¬ 
trollability  unless  the  control  input  is  distributed  along  the 
entire  interval  in  the  physical  domain.  In  the  presence  of 
inertia,  crucial  difference  occurs  between  the  cases  of  one 
or  more  relaxation  modes:  For  a  single  relaxation  mode 
(i.e.,  a  Maxwell  fluid),  exact  controllability  holds  provided 
the  time  interval  satisfies  certain  inequality;  for  multiple 
relaxation  modes,  exact  controllability  holds  under  mod¬ 
ified  regularity  assumptions.  Another  piece  of  work  by 
Renardy90  was  focused  on  the  homogeneous  shear  flow  of 
viscoelastic  fluids  with  several  different  constitutive  mod¬ 
els.  For  those  equations,  the  state  of  the  system  consists  of 
viscoelastic  stresses  whereas  the  shear  rate  is  regarded  as 
a  control  input.  For  the  upper  convected  Maxwell  (UCM) 
model,  it  was  revealed  in  Ref.  [90]  that  the  reachable  set, 
i.e.,  the  states  in  stress  space  which  are  accessible  from  a 
given  initial  condition,  is  specified  by  a  positive  definite¬ 
ness  inequality  of  the  stress  tensor.  Very  recently  the  result 
was  extended  to  the  control  of  nonhomogeneous  shear  flow 
of  an  upper  convected  Maxwell  fluid91  where  the  state  of 
the  system  and  the  available  control  are  the  same  as  those  in 
Ref.  [89].  In  our  very  recent  studies121  we  started  to  apply 
nonlinear  geometric  control  theory  to  examine  the  control¬ 
lability  of  the  upper  convected  Maxwell  model  under  vari¬ 
ous  homogeneous  flows.  In  Ref.  [121]  we  adopted  a  local 
version  of  the  concept,  namely  weak  controllability.  This 
definition  has  been  widely  used  in  nonlinear  control  the¬ 
ory  as  well  as  in  engineering  applications.  In  many  real  life 
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applications  of  complicated  nonlinear  systems,  engineers 
prefer  to  reach  a  distance  target  by  scheduling  a  sequence 
of  local  movements  for  the  reasons  of  model  uncertainties, 
system  perturbations,  and  sensor  noise. 

In  the  following  we  present  on  overall  basic  defini¬ 
tion  and  description  of  weak  controllability  adopted  from 
Ref.  [55], 

Let 

m 

*  =  /(•*)  + 1] &(XK  (44) 

;=i 

be  a  general  nonlinear  control  system,  where  x  is  the  state 
variable,  and  m,  e  SK ,  i  =  1, . . . ,  m,  are  the  control  variables. 
Let  M  denote  the  manifold  of  state  variables.  A  point  x, 
in  M  is  reachable  from  a  point  x0  in  M  if  there  exist 
piecewise  continuous  input  functions,  tq  =  a, -(f),  such  that 
the  trajectory,  x(t),  of  (44)  with  initial  state  x0  reaches  x, 
in  finite  time.  The  global  reachability  is  usually  hard  to 
prove  for  nonlinear  control  systems.  A  feasible  solution 
is  to  seek  weak  controllability.  The  system  (44)  is  weakly 
controllable  within  some  open  subset  S  C  M  if  for  each 
point  x0  e  S,  there  is  an  open  neighborhood  U0  of  x0  so 
that  the  set  of  points  reachable  from  x0  along  trajectories 
inside  U0  contains  at  least  an  open  subset  of  M. 

A  powerful  sufficient  condition  on  the  weak  controlla¬ 
bility  of  a  nonlinear  control  system  is  the  contollability 
rank  condition  (CRC).55  Namely,  a  control  system  of  the 
form  (44)  is  weakly  controllable  on  an  open  set  S  if  it 
satisfies  the  controllability  rank  condition  on  S,  i.e., 

dim(  Acg(x))  =  n  (45) 

for  all  x  e  .S',  where  n  is  the  dimension  of  the  manifold  M, 
and 


Acg(x)  =  span{X(x)  |  Xis  a  vector  field 

in  the  control  Lie  algebra}. 


From  a  theoretical  viewpoint,  Lie  brackets  of  vector 
fields  provide  an  efficient  tool  to  prove  weak  controllability. 

The  study  in  Ref.  [121]  can  be  extended  to  the  nematic 
liquid  crystal  polymers  using  the  Doi-Hess  tensor  theory. 
As  a  toy  model,  let  us  consider  the  two-dimensional  Doi 
model  with  imposed  extensional  flow. 

For  a  homogeneous  extensional  flow  with  rate  y  (r),  the 
velocity  is 

v  =  j»(|.-0  (46) 

The  rate-of-strain  tensor  becomes 


D=I[Vv+(Vv)J']=^ 


1 

0 


(47) 


The  Doi  tensor  model  for  nematic  liquid  crystalline  poly¬ 
mer  under  extensional  flow  is  given  by 

d 

—  Q  =  flQ-Qfl  +  a[DQ  +  QD] 
dt 

+  oD  —  2aD  :  Q(Q  + 1/2)  —  6DrF(Q)  (48) 


where 


Gn  Qi 


0i 


-0i 


is  the  orientation  tensor, 


1  rp 

Cl  —  -[Vv  —  (Vv)  ]  is  the  vorticity  tensor, 

F(Q)=  (i-  y)Q-NQ2  +  AQ:Q(Q+* 


(49) 


a  is  a  dimensionless  parameter  depending  on  the  molecular 
aspect  ratio,  A  is  a  dimensionless  concentration  of  nematic 
polymers,  Dr  is  an  averged  rotary  diffusivity  or  relaxation 
rate.  After  some  rescaling  and  nondimensionalization,  the 
nematodynamic  model  (48)  can  be  written  as 


where 


dx 

dt 


=  f(x)  +  g(x)u 


1 

_ 1 

"  011  " 

|_*2_ 

.  012  . 

fix)  =  -6 


1  — —  +  2A(.Xj  +  x2) 


—2axxx2 


u  =  Pe(t) 


(50) 


We  now  regard  Pe(t)  as  a  control  input  and  consider  the 
problem  (50)  with  initial  state  Q0  and  final  state  Q  .  The 
system  is  controllable  if  for  any  choices  of  Q0  and  Q, . 
there  exists  a  control  Pe(t)  on  the  interval  (0,  T). 

The  Lie  bracket  is  calculated  to  be 


6  a 

\ 2x 2 

N 

1--- 2N(x2{  +  x2) 

II 

fSo 

IN  ,  1 

+  -  —  —  +  A  (3x[  +  x2)  | 

12flX!X2 

i  i 

<N  (N 

+ 

(N  r-H 

(N 

1 

^  |<N 
+ 

1 _ 1 

Finally,  we  have 

det  [g,  [/,  gJJ  =  12a2X|X2  (51) 


The  determinant  does  not  equal  zero  when  x,  ^  0  and 
x2  ^  0.  So  the  Doi  model  (50)  under  extensional  flow  sat¬ 
isfies  CRC  and  is  weakly  controllable  when  Qn  ^  0  and 
Qn  ^  0.  In  terms  of  geometric  words,  the  system  (50)  is 
weakly  controllable  at  all  points  in  R2\{5’1  US2}  where 

R2  =  { (xj , x2)  |  X|,x2  e  R},  .5,  =  {(x,,x2)  |  xt  =0}, 

5’2  =  {(x1,x2)  |  x2  =  0} 

Extending  the  above  preliminary  study  to  3-D  orienta¬ 
tion  tensor  theory  and  other  flow  geometries  (e.g.,  shear 
flow)  and  finding  reachable  set  remain  to  be  explored. 
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The  observability  of  nematic  liquid  crystal  polymers  is 
still  open.  Roughly  speaking,  by  observability  we  mean 
the  ability  to  determine  the  orientation  distribution  or  ori¬ 
entation  tensor  of  nematic  polymers  from  the  time  his¬ 
tory  of  the  observations.  The  observations  can  vary,  for 
example,  the  birefringence  measurements  from  light  scat¬ 
tering  experiments.  Recently  Krener65  studied  the  observ¬ 
ability  of  one  and  two-point  two-dimensional  vortex  flow 
from  one  or  two  Eulerian  and  Lagrangian  observations. 
In  his  work,  an  Eulerian  observation  is  a  measurement  of 
the  velocity  of  the  flow  at  a  fixed  point  in  the  domain 
of  the  flow  whereas  a  Lagrangian  observation  is  a  mea¬ 
surement  of  the  position  of  a  particle  moving  with  the 
fluid.  By  applying  the  observability  and  strong  observ¬ 
ability  rank  conditions  and  calculating  them  for  different 
vortex  configurations  and  observations,  it  was  found  that 
vortex  flows  with  Lagrangian  observations  tend  to  be  more 
observable  than  the  same  flows  with  Eulerian  observations. 
Extending  this  work  on  vortex  flows  to  complex  flows 
including  nematic  liquid  crystal  polymers  is  practically 
important. 

7.  OTHER  DIRECTIONS 

Other  current  and  future  research  directions  of  nematic 
liquid  crystals  include  theoretical  and  numerical  studies 
of  biaxial  liquid  crystals  whose  molecules  are  ellipsoidal, 
brick-shaped,  V-shaped,  bow-shaped,  boomerang-shaped 
or  banana  shaped.76  Recall  that  in  a  nematic  phase  the 
molecules  tend  to  align  along  the  director  n.  A  biax¬ 
ial  nematic  phase  may  result  due  to  the  further  break¬ 
ing  of  the  rotational  symmetry  of  the  system  around  the 
director  n.  The  existence  of  the  biaxial  nematic  phase 
was  originally  predicted  theoretically45  and  later  observed 
experimentally.78 

The  equilibrium  phase  diagram  of  a  class  of  biaxial 
nematic  liquid  crystals  was  given  in  Refs.  [5, 99]  using 
mean  field  theories  with  a  generalized  Straley’s  interac¬ 
tion  potential.  Very  recently  shear  induced  mesostructures 
(either  steady  or  time-dependent)  in  biaxial  liquid  crystals 
were  reported  in  Ref.  [97]  using  a  hydrodynamical  kinetic 
theory  where  the  governing  Smoluchowski  equation  was 
solved  numerically  with  the  Galerkin  method.  In  Ref.  [97] 
the  simple  plane  shear  flows  were  used.  Extending  this  to 
more  complicated  flows  remains  to  be  explored. 

8.  CONCLUDING  REMARKS 

In  this  review  we  have  briefly  summarized  recent  theo¬ 
retical  and  numerical  studies  on  liquid  crystal  polymers. 
Some  unexplored  problems  and  areas  which  need  future 
attention  have  been  identified.  Due  to  its  complexity  and 
importance,  this  field  will  inspire  and  involve  more  scien¬ 
tific  activities. 
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